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ABSTRACT 

Using high-quality Hubble Space Telescope observations, we construct the near infra-red 
(NIR) to far ultra-violet (FUV) spectral energy distribution (SED) of PSR B0656-hl4. The 
SED is non-monotonic. Fitting it with a simple combination of a Rayleigh- Jeans spectrum (UV) 
and non-thermal power-law (optical/NIR) leaves significant residuals, strongly hinting at one 
or more spectral features. We consider various models (combination of continuum components, 
and absorption/emission lines) with possible interpretations, and place them in the context of the 
broader spectral energy distribution. Surprisingly, the extrapolation of the best-fit X-ray spectral 
model roughly match the NIR-FUV data, and the power-law component is also consistent with 
the 7-ray fluxes. We compare the multiwavelength SED of B0656-f 14 with those of other optical. 
X-ray and 7-ray detected pulsars, and notice that a simple power-law spectrum crudely accounts 
for most of the non-thermal emission. 



Subject headings: pulsars: individual (PSR B0656+14) 
1. Introduction 

Pulsars are responsible for the production of 
high-energy emission observable across the EM 
spectrum from the radio to 7-rays. They are nat- 
ural laboratories for investigating particle acceler- 
ation and radiation processes in the high-energy, 
high-field, high-gravity regime. Outside the radio, 
the non-thermal high energy emission is believed 
to be produced by incoherent processes such as 
synchrotron and curvature radiation. With cur- 
rent instrumentation, this radiation is most eas- 
ily detected in 7-rays (where most of the energy 
is released) and X-rays (although the latter is 
often "contaminated" by thermal emission from 
the neutron star surface). It is much more chal- 
lenging to detect the non-thermal emission in the 
optical / near-infrared (NIR) . 

There are few pulsars for w hich mul t icol or 
UV/optical/IR photometry exists (Mign ani 200^ 
Kargaltsev fc Pavlov! [2007t iMignani et al.. ,20071) . 
These are moderately old pulsars B1929-I-10 and 
B0950-f08; r ~ 3Myr and 17.5 Myr, respectively). 



middle-aged pulsars Geminga, B0656+14, and 
B1055-52 (~ 300, 100, and 500 kyr respectively), 
the younger Vela pulsar (~ 10 kyr), and the 
very young pulsars Crab (950 yrs) and B0540— 69 
(~ 1.7 kyr). The optical and X-ray spectra of the 
Vela, Crab and B0540-69, look "boring" (feature- 
less power-laws), while the spectra of the middle- 
aged pulsars are more complex and exhibit both 
thermal (emitted from the NS surface) and non- 
thermal (magnetosphe ric) components in X-rays 



and optical-UV (e.g., Romani et al.l lioosl ) . No 



spectral feature has yet been observe d in the IR- 
UV sp ectra of pulsar^, although iZharikov et al 



(|2007l )'s data hinted at possible feature(s) for PSR 
B0656-I-14. Here we will present detailed analysis 
of the IR-UV spectrum of B0656-f 14. 

PSR B0656+14 (P = 385 ms; t = 1.1 x 
lO'^yr), located at the distance of 288 ± 30 pc 



teignami et at I l ll996fl claimed a feature in the optical- 
UV spectrum of Geminga (see also iMignani et al.l Il998l ) , 
but this was not supp orted by subsequent observations 



but this was not supp 
jKargaltsev et al.ll200g) 
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(jBrisken et al.l 120031 ) , is the brightest of middle- 
aged pulsars, in both X-rays and optical. Its 



polarized. Finally, the low-resolution spectrum 
and pulsations in the near- UV (NUV) have bee n 



by Anderson et al. (1993): 


Possenti et al. (1996): 


Greiveldinger et al. ( 1996 


); Marshall & Schulz 


(2002 


): Zavlin & Pavlov ( 


2004: iDe Luca et al. 


(2005 


). The X-ray spectrum (0.1-6 keV), ob- 



2002), can be fitted with a model which consists 
of thermal soft (TS) component with a blackbody 
(BB) temperature of ~ 0.82 MK, emitted from a 
large part of the NS surface (i? ~ 7.3 km); thermal 
hard (TH) component with the BB temperature 
of ~ 1.7 MK, apparently emitted from smaller 
{R ~ 0.5 km) hotter areas (perhaps hot polar 
caps); and non-thermal power-law (PL) compo- 
nent, with the photon index F ~ 1.5, possibly 
emitted from the pulsar magnetosphere. The X- 
ray radiation is pulsed, with a pulsed fraction of 
~ 13% at lower energies (< 0.7 keV), where the 
TS component dominates, and ^ 57% in the 2 — 4 
keV band, where the PL component is dominant. 
The optical counterpart of P SR B0656-H4 
was discovered by ICaraveo et al 



(X. 



25) 



(|1994| ). Because a near by galaxy contaminate s 
the ground-based images (jKoptsevich et al ] l200ll) . 
the best near-IR/optical/near-UV (A = 0.2- 
2 /xm) data were obtained with the Hubble 
Space Teles cope (HST). The p u lsar w as imaged 



with FOG (IPavlov et al.l Il996l IiQQTI). WFPC2 



( Mignani et al.ll2000D . and NICMOS (|Koptsevich et al, 
20011 ) . These observations suggested that the pul- 
sar has a non-thermal spectrum from NIR through 
the optical, with somewhat uncertain spectral in- 
dex, ~ -0.4 ± 0.5 {F^ oc u""). The large un- 
certainty in the slope reflects a substantial scatter 
among the existing photometric points. This scat- 
ter could possibly be attributed to the presence of 
broad spectral features or could be just the result 
of poor cross-calibration between different instru- 
ments that have been used to cover a broad range 
of wavelengths. 



In a ddition to photometric observations . iKern et al 
(12003) observed the optical pulse profile and po- 
larization with the Palomar 5 m telescope and 
found that the emission in 400— 600 nm band 
is highly pulsed (upper limit on the un-pulsed 
flux is 16%). The pulse profile is double-peaked, 
with a bridge of emission between the two peaks; 
emission from the bridge was found to be highly 



200i). 



measured with HST STIS (|Shibanov et al 
They found that NUV pulse profile is also double- 
peaked and similar to the optical one. The high 
pulsed fraction, 70 ± 12 %, and a flat spectrum 
(tti/.NUV = 0.35 ± 0.5) suggest that NUV emis- 
sion is mostly non-thermal. Since the soft X-ray 
spectrum is thermal and optical spectrum is non- 
thermal, a transition between the thermal and 
non-thermal spectra is expected to occur some- 



where in between. iKargaltsev fc Pavlov! ([20071) 



show that the FUV flux is consistent with a low 
pulsed fraction thermal component and a highly 
pulsed non-thermal component. 

To better constrain the shape of the spectrum 
and evaluate contributions of the thermal and non- 
thermal components, we observed B0656 using in- 
struments aboard HST (ACS, STIS and COS). To- 
gether with re-analysis of archival NICMOS obser- 
vations, we present the most complete infrared to 
ultraviolet spectral energy distribution measured 
for any pulsar. The data analysis of the obser- 
vations, including the archival ones, is described 
in Section 2 and Appendix A. In Section 3, we 
present the results, both the photometry and spec- 
tra from the individual observations, and the com- 
bined spectral energy distribution. In Section 4 we 
discuss the interpretation of the measured spec- 
tral shape, its context in the wider spectral energy 
distribution, and comparisons with the spectra of 
other pulsars. 

2. Observations and analysis 

In Table [T] we list the observations analyzed as 
part of this work. Of these, the NICMOS and 
STIS NUV/FUV observations were reported be- 
fore (iKoptsevich et aLlbOOlH Shiba nov et al.| [200l 
Kargaltsev fc Pavlovl 120071 ). but we reanalyze the 
data from the archive with the latest calibration 
applied. The rest of the data appear here for the 
first time. 

Table [2] summarizes previously published opti- 
cal/UV observations of B0656-I-14, which we will 
not re-analyze in this paper. We will include the 
results of these observations with those produced 
in this work, for comparison. 

Below, we briefly describe the fluxes and spec- 
tra derived for each data-set, in order of increasing 
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Table 1 

Log of HST observations for analysis 





Frogrctm ID 


instrument / 
optica,! elGnient 


Wcivclcngtn 


Exposure 
tiriiG (s) 


1998-03-18 


7836 


NICMOS/FllOW 


11285 


2544 




7836 


NICMOS/F160W 


16056 


5089 




7836 


NICMOS/F187W 


18718 


7633 


2001-09-01 


9156 


STIS/prism 


1750-3000 


6791 


2001-11-16 


9156 


STIS/prism 


1750-3000 


12761 


2004-01-20 


9797 


STIS/G140L 


1150-1700 


4950 


2005-12-03 


10600 


ACS/FR647M 


6820 


810 




10600 


ACS/FR647M 


7393 


810 




10600 


ACS/FR914M 


7956 


747 




10600 


ACS/FR914M 


8751 


1331 




10600 


ACS/FR914M 


9491 


2842 


2005-12-08 


10600 


ACS/FR459M 


4029 


1400 




10600 


ACS/FR459M 


4401 


1020 




10600 


ACS/FR459M 


4780 


830 




10600 


ACS/FR459M 


5165 


885 




10600 


ACS/FR647M 


5566 


1020 




10600 


ACS/FR647M 


6244 


1020 


2010-02-17 


11629 


COS/G140L 


1100-1700 


10927 


2010-02-20 


11629 


COS/G140L 


1100-1700 


10927 



Note. — Wavelengths are either the sensitive range of a spectroscopic observa- 
tion, or the pivot wavelength of a photometric observation. 



Table 2: List of previously published observations and results. 



Instrument Filter Central Frequency Flux F^, Reference 

(10" Hz) (/xJy) 



HST/WFPC2 


F555W 


5.48 


0.39(2) 


1 


HST/FOC 


F555W 


5.62 


0.39(2) 


2 


HST/FOC 


F430W 


7.28 


0.26(3) 


3 


HST/FOC 


F342W 


8.81 


0.31(3) 


3 


HST/FOC 


F195W 


12.7 


0.35(4) 


3 


BTA/CCD 


I 


3.80 


0.60(6)" 


4 


BTA/CCD 


R 


4.61 


0.39(3) 


4 


BTA/CCD 


V 


5.46 


0.37(4) 


4 


BTA/CCD 


B 


6.84 


0.32(4) 


4 


Subaru/Suprime 


I 


3.80 


0.37(4) 


5 


Subaru/Suprime 


R 


4.61 


0.42(3) 


5 


Subaru/Suprime 


B 


6.84 


0.31(2) 


5 


VLT/F0RS2 


4300-9600 A 






6 



Note: references (l ) lMignani et~al1 (Il997h (2) |Pavlov et al.l (Il99e U3) IPavlov et al.l (Il997l) (4) 
Koptsevich et al. I (|200lh (5) IShibanov et al.l (12006^) (6) IZharikov et all (|2007h . Numbers in parentheses 



indicate the uncertainty in the final digit. 
The authors note that this measurement was contaminated by a nearby red galaxy. 
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wavelength. Full details of the data reduction pro- 
cedures used are given in the Appendix. We en- 
courage interested readers to refer to the Appendix 
for a thorough account of the specific procedures 
required for the analysis of the diverse data, in par- 
ticular grating spectroscopy (with a slit), slit-less 
prism spectroscopy, very faint infrared photome- 
try (requiring custom NICMOS darks), ramp filter 
photometry and the evaluation of calibration per- 
formance in the IR-optical regime. 

2.1. STIS FUV grating 

The measured flux values are given in Ta- 
ble [21 while the spectrum is shown in Fig- 
ure [B The total flux in the 1153-1700 A 
range (AA = 547 A), can be estimated as F ~ 
AA (E^(^A)^AAO (E. AA,)"' ^ (4.28 ± 0.28) x 
10~^^ erg s~^ cm~^, corresponding to the lumi- 
nosity Lfuv = ^tt(PF ^ (4.24 ± 0.28) x 10^8^2^^ 
erg s~^. 

For a typical NS radius R — 13 km and distance 
d — 288 pc, the inferred brightness temperatures 
are 0.615 ±0.026, 0.663±0.028, 0.715±0.031, and 
0.832 ± 0.036 MK, for E{B - V) = 0.01, 0.02, 
0.03, and 0.05, respectively; the corresponding 
values are 0.663, 0.642, 0.626, and 0.611, for 4 
degrees of freedom (dof). An example of best- 
fit blackbody spectrum is shown in Figure [U for 
a plausible E{B — V) = 0.02. The correspond- 
ing unabsorbed bolometric luminosity can be es- 
timated as Lboi = 1-2 X Kfi^T^Rl^ erg s'^ for 
instance, iboi = (2.3 ± 0.4) x 10^^ erg s"^ for 
E{B ~V) = 0.02, i? = 13 km and d = 288 pc. 

2.2. COS spectroscopy 

The measured spectrum is shown in Figure 
[2] while the spectral fluxes are listed in Table 
H An absorbed PL, Fx oc A""\ with best-fit 
a A = 3.0 ± 0.2 is plotted. The curve follows the 
overall shape of the spectrum well, but there is 
substantial scatter resulting in = 2.1. The ma- 
jority of this scatter is in a few outlier points which 
not appear to form spectral feature(s). We cannot 
rule out narrow lines, but such lines are not ex- 
pected for the high temperature and pressure of 
a NS atmosphere or magnetosphere. We find no 
obvious reason for the additional scatter, and spec- 
ulate that it may be due to the imperfections of 
the calibration or structure in the detector back- 



ground. 

2.3. STIS NUV prism 

The measured flux values are given in Table [5l 
while the spectrum is shown in Figure[3l The total 
flux in the 1790-2950 A range (AA = 1160 A), can 
be estimated as (2.63±0.38) x 10"^^ erg cm^^ 
corresponding to the luminosity Lnuv = 4:Trd^F = 
(2.61 ± 0.38) x lO^^d^gg erg s'^ The obtained 
average flux is slightly larger, but compatible with 
the flux measured by Shibanov et al. (2005) in 
somewhat wider wavelength interval (1700 — 3400 

A). 

We then fit the NUV spectrum with the ab- 
sorbed power- law model, F\ = F2500 (A/2500 
A)"A X io-oAAix)E{B-v)_ plausible values 

E{B -V)^ 0.01, 0.02, 0.03 and 0.05, we find the 
power-law indices ax = -2.92±0.41, -3.01±0.41, 
-3.09 ±0.41 and -3.24 ±0.41, and the normaliza- 
tion F2500 = 1-85 ± 0.08, 1.97 ± 0.09, 2.11 ± 0.10 
and 2.41 ± 0.11 x lO^^^ erg cm'^ s'^ A^^, re- 
spectively; the corresponding xt values are 1.27, 
1.38, 1.51, and 1.87, for 5 dof. The obtained spec- 
tral slopes are systematically smaller than those 
derived from our fits to the FUV spectrum. The 
flattening of the spectrum at longer wavelengths 
can be attributed to a larger relative contribu- 
tion of the non-thermal component. The measured 
slope of the NUV spectrum for E{B -V) ^ 0.03, 
~ 1.1 ± 0.4 (tty = —ax — 2), is marginally com- 
patible with the slope a^, = 0.35 ± 0.5 found by 
Shibanov et al. (2005). The difference can be at- 
tributed to the different choice of the spectral bins 
and narrower wavelength range that we used. 

2.4. ACS/WFC Ramp Filter Photometry 

The photometric fluxes for three ramp filters 
and a total of 11 central wavelength positions are 
listed in Tableland shown in Figure |4l Given the 
non-standard nature of the ramp filter photome- 
try, we performed a specific check on our calibra- 
tion in Section IA.6I We find the pipeline calibra- 
tion satisfactory. 

2.5. NICMOS photometry 

We find F^(FllOW) = 0.385 ± OMOfi-fy; 
F^(F160W) = 0.551 ± 0.030 AiJy; Fj.(F187W) = 
0.660 ± 0.025 /iJy (uncertainties are 1-cr and sta- 
tistical only). Due to the non-standard processing 
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I I I I I I I I I I 

1100 1200 1300 1400 1500 1600 1700 

Wavelength, A 



Fig. 1. — STIS FUV spectrum of B0656. The solid line shows an absorbed black-body spectrum with 
T = 0.663 MK and E{B -V)= 0.02. 
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1200 1400 1600 

A (A) 

Fig. 2. — COS-FUV spectrum of PSR B0656+14. The gaps are due to the discarded data strongly 
affected by geocoronal emission hues (see Section 12.21) . An absorbed PL fit is plotted, with Fx = 
A-"^ X lO-oAA(x)EiB-v) ^ ^ 3 ± 0.2, E{B -V)^ 0.03. 



6 



steps and calibration involved, we conservatively 
estimate the systematic uncertainty to be about 
10% in each case (see Section IA.6j) . Our flux 
values are o nly marginally consisten t with those 
reported bv iKoptsevich et aL ( 2001 ). We believe 
that by using the most up-to-date calibration, our 
values should be the most accurate. We verify our 
calibration in Section lA. 61 

3. Results and discussion 
3.1. IR-FUV spectrum 

The final spectrum shown in Figure S] includes 
both the fluxes/spectra calculated in this work and 
those from the literature (see Tabled]). This is, to 
our knowledge, the most complete IR-UV spec- 
trum of a non-recycled pulsar. 

3.1.1. Qualitative spectral shape 

Clearly, the spectrum of B0656 shown in Figure 
IHdoes not fit to a simple absorbed PL model. The 
spectrum appears to be more complex, with some 
structure. Therefore we begin by making some 
qualitative statements. 

The upturn of the spectrum in the NUV range 
is consistent with a contribution by a thermal, 
Rayleigh- Jeans (R-J; oc v^) component that 
is most prominent in FUV and also consistent 
with the extrapolation of t he X-ray TS compo- 
nent to the optical (see also iKargaltsev fc Pavlov 
20071 ). However, even with the thermal compo- 



nent added, any simple PL fit to the optical spec- 
trum of B0656 leaves large residuals. In particu- 
lar, there is a deep minimum in the spectral fiux 
around v = ?> x 10^"* Hz (A = 1/im), with steep 
rises to either side (see Figure H]). 

The complex spectral shape can be interpreted 
in several ways. For instance, the steep rise at the 
lowest frequencies could be attributed to a cool, 
T ~ 1000 K, thermal (disk) component. Alter- 
natively, the spectrum is very suggestive of spec- 
tral features^ either two absorption lines or a very 
broad emission line. We have fitted models cor- 
responding to these three scenarios with the best- 
fit parameters listed in Table [T] The correspond- 
ing models and residuals are shown in Figure [5] 
We note that none of the fits is statistically good 
(x^ > 1.5), but these fits do not show system- 
atic structure in the residuals. Of the three. 



the absorption model, with lines at (2.91±0.03) 
and (8.05±0.09) x lO^Hz (A = 1.030±0.011 and 
0.372±0.004/im), is statistically preferred. 

3.1.2. Possible interpretations 

Several possible interpretations for spectral fea- 
tures can be considered. As an example, we illus- 
trate the apparently straightforward case of cy- 
clotron emission (or absorption). Given the cur- 
rent uncertainties of the optical flux, we could 
speculate that we are observing cyclotron lines 
from a cool plasma. The NS magnetosphere has 
a dipole magnetic field B cx [RNs/f)^ with sur- 
face strength iJdipoie 5 x 10^^ G . An electron 
cyclotron line centred at A w 1 /im would require 
magnetic field i3 '--^ 2 x 10^ G, which would place 
the line forming region at d « ZOORns from the 
NS surface but still well within the light cylinder 
radius (i?LC = 1.15 x IO^Ens)- 

A deep dip at 1/im (i^ = 3 x 10^"* Hz; Figure 
|4|) could signify a strong absorption line, and so 
suggest the hypothesis of an absorbing cyclotron- 
resonance layer which could b e composed of e"'"/e~ 
(Raiagopal & Romanil Il997t IWang et al.l Il998t 
Jludcrman 2003). Such a layer can significantly 
reprocess/change thermal radiation from the hot 
NS surface. If confirmed by more sensitive ob- 
servations of B0656, quantitative measurements 
(centroid, line width, phase dependence) will open 
a new window to probe the structure of the pulsar 
magnetosphere. 

Alternatively, there may be a cold thermal com- 
ponent in the infra-red (a steepening rise towards 
the red) which could be indi cative of a pas sive 
post-supernova fall-back disk (|Chevalier Il989l) for 
which we estimate a temperature of about 500- 
1500 K. This would suggest a link to magnetars, 
where s uch fall-back disks have already been sug- 
gested (|Wang et al.l l2006l) . and support recent 



models that predict the persistence of fall-back 
disks to much later epoch s than recently thought 
(jD'Angelo fc SpruitJl201lh . 



3.2. Broad-band multi-v^ravelength spec- 
trum 

Let us consider the NIR-UV spectrum of B0656 
in a broader, multi-wavelength context. B0656 
has been extensively studied in X-r ays (e.g., 
Pavlov et al.l 120021 : iDe Luca et al.l l2005l ) , and has 
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also been detected as a 7-ray source by Fer mi LAT 
with high significance (jAbdo et ah I l2010bh . 



We downloaded Chandra ObsID 2800 (25 ks 
exposure with ACIS in continuous clocking [CC] 
mode; see also Pavlov et al.l 1200*2 ) and also the 
spectrum of ObsID 130 (38 ks exposure with 
the HRC a nd low-energy transmiss ion grating 
[LETG]: see lMarshah fc Schuljl2002( ). The data 
were reduced using standard procedures. We con- 
vert the counts spectra to fluxes by dividing the 
count rate in each bin by the calibrated effective 
area in that bin (sometimes called ARF), after 
convolving with the response function (sometimes 
called RMF). Although such 'fluxing' can produce 
an incorrect spectral shape at sharp detector sensi- 
tivity edges and spectral featurefjj, the deviations 
are confined to relatively narrow energy ranges, 
and the spectrum obtained in this way is inde- 
pendent of any pre-supposed spectral model. We 
also use d the fluxes from th e second Fermi catalog 
(2FGL; iThe Fermi-LAT Collabora tion 2011 ) . 

In Figure IHl we show these data along with the 
NIR-UV fluxes. Note, however, that the extinc- 
tion (and optical/UV reddening) is rather uncer- 
tain, and has a particularly strong effect at the 
softest X-ray energies {E ~ 100 eV). 

We have plotted in Fi gure [SI the best- fit X-ray 
model from lPaylov et al.l (|2002') . fitted to the same 
Chandra dat^f[ The extrapolation of the X-ray 
model, the sum of two black-bodies and a power- 
law, is shown in the inset of Figure |6l Surprisingly, 
both components fall close to the UV-optical spec- 
trum (Section l3.ip . Extrapolating the X-ray fits 
into the NIR-UV exacerbates the uncertainties in 
the model flux: a small uncertainty in any given 
parameter translates into a large uncertainty in 
the predicted fiux in the NIR-UV regime. It is nat- 
ural to assume that the Rayleigh- Jeans spectral 
component is a continuation of the cooler X-ray 



^X-ray detectors are prone to 'redistribution', the possibility 
that a photon produces a signal corresponding to a different 
energy. 

^De Luca et al.l ||2005| ') fitted the same model to XMM- 
Newton data, and found fit parameters formally incon- 
sistent w ith the ones obtained from the Chandra data; 
whereas l|Zavlin fc Pavlovl |200A found consistent param- 
eters. We note that even within the XMM-JVewton data 
there are obvious discrepancies between the spectra from 
two differe nt instrument modes for E < 0.5 keV (see Fig- 
ure 2 from lOe Luca et al.ll20d5l ) and therefore the best-fit 
uncertainties have likely been underestimated. 



black-body; for a temperature of T = 0.82 MK, its 
radius would be i? ~ 10 km (as observed at infin- 
ity) to fit the FUV, with a 15% uncertainty on T 
for the reddening range 0.01< E{B - V) <0.05. 
The power-law component in the optical could, in 
principle, be a continuation of non-thermal power- 
law emission seen in X-rays (but would not be suf- 
ficient to describe the data without the features or 
a cool disk; see Section [3. 1.1|) . 

The extrapolation of the best-fit PL also hap- 
pens to fall close to the Fermi 7-rays fiuxes, hint- 
ing that a single very broad PL with photon in- 
dex Fmw = 1.50 ±0.05 could in principle describe 
the entire non-thermal spectrum from optical to 
7-rays. In order to match the GeV fiuxes, the PL 
should of course be modified by a high -energy cut 



off, no rmally seen for 7-ray pulsars (l Abdo et al 
2010bl ). The high pulse d fraction (>50%) in the 



Fermi LAT hght- curve (lAbdo et al.l l2010bl ) and 



also in the optical ( Kern et al.ll2003ir confirms that 
both are magnetospheric. 

3.3. Comparison to other pulsars 

Less than 1% of radio pulsars have been de- 
tected in the NIR-optical-UV spectral window to 
date. Of those all have been detected in X-rays 
but only some in 7-rays (Table We note that 
this sample includes members that are young and 
energetic, middle-age pulsars such as B0656+14, 
and one millisecond (recycled) pulsar (MSPs) . No 
isolated (i.e., non-accreting) pulsar has yet been 
discovered in the IR-UV range without having pre- 
viously been detected in the radio and X-rays. 

We retrieved the IR/optical/UV fluxes and X- 
ray spectra of the pulsars in TablelH using our pre- 
vious work (e.g., Durant et al., 2011: |Pavlov et aT 



12001: Zavlin fc Pavlov.2004.) . X-ray spectral prod- 
ucts downloaded from the Xassist projecl|f|, or, for 
the case of B1951-I-32, by re-processing data from 
the archives. Note that the l ow-significance op- 



tical detection of B1951-I-32 ( Butler et al.l 12002 : 



Moon et al.l 120041 ) has not been confirmed. The 
GeV fluxes were taken from the 2F LG catalog 
()The Fermi-LAT Collaboration! [2QTTh . 

In Figure |7| we show the IR-GeV spectra of 
the pulsars in Tab le [U except for the Crab (see 
Kuiper et all l200l[) . Fi gure El shows the same 



' ^http: //xassist .pha. jhu. edu/zope/xassist I 
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data, but in terms of 'spectral efficiency', 77^ = 
Att(PvFv/E. We have sketched PL hnes which 
connect the 7-ray fluxes with the optical for each 
spectrum and find that in each case the line also 
passes through the high-energy, non-thermal part 
of the X-ray spectrum. This suggests that each 
pulsar's non-thermal spectrum can be empirically 
described by a single PL stretching across many 
orders of magnitude in frequency, w ith a break 
or cut-o ff at GeV energies. Note, that lAbdo et al 



(|2010a[) have already considered whether the 7-ray 
and X-ray non-thermal spectra may be consistent, 
and conclude that they are not (see their Figure 
5); however, there is no discussion there of the un- 
certainties in the fitted parameters or of the effect 
of the choice of model function (e.g., atmosphere 
models for X-rays and broken power-laws in 7- 
rays). In Table [8] we give the appropriate spectral 
indices of the lines we have sketched (these are not 
true fits to the data). The range of the photon in- 
dex, 1.2 < Fmw < 1-65, is remarkably narrow. We 
note that the scales in the Figure are logarithmic, 
and there may well be additional structure on top 
of the speculative underlying PL component (as 
detected in B0656, see Figure |6|), but surprisingly 
it appears that a remarkably simple description 
can be applied for a number of pulsars. We note 
that such a simple description is not universal, as 
it is not appropriate for very young pulsars such 
as the CrahH. 

Recentlv iDanilenko et al. ( 2011 ) presented ev- 
idence for steeply-rising infra-red excesses in the 
Geminga and Vela pulsars (the former only a 
marginal detection) . It may be that this is a com- 
mon feature of middle-aged pulsars. For Geminga, 
which has a similar spectrum, luminosity and dis- 
tance to PSR B0656-I-14, they conclude that the 
IR excess (if it is real) can be successfully modeled 
with a cool irradiated disk truncated at an inner 
radius comparable to the light cylinder Rlc = 
O.O16i?0. In the case of PSR B0656-I-14, we do 
not have good enough data to attempt to model 
disk- like emission, and Spitzer observations cannot 
resolve the pulsar from the brighter nearby galaxy. 
Ground-based adaptive optics in the 2-4 /im range 
and/or better signal-to-noise in the current infra- 



^The recent high-energy spectrum of the Crab is shown in 
the recent poster confluence.slac.staiiford.edu/ 
download/attacliments/102860834/ 
KUIPER_PSR_S2 . N19_PDSTER . pdf 



red range (using the newer Wide Field Camera 3 
on HST) should clarify the situation. 

4. Summary 

The IR-UV spectrum of PSR B0656-hl4 is non- 
monotonic, and requires a combination of emission 
or absorption lines and/or continuum components. 
The putative features might be explained as cy- 
clotron lines, hinting at an outer-magnetospheric 
electron belt. There is also the possibility of a 
cold (T ^ 1000 K) thermal component, such as a 
passive, dusty circumstellar disk. Further observa- 
tions in the IR are required to confirm and probe 
the potential disk component. 

The multi- wavelength spectrum of PSR B0656-I-14 
suggests that a single, very broad PL component 
can describe the non-thermal emission surprisingly 
well. Interestingly, there may be evidence for such 
a component in other pulsars that have both 7-ray 
and IR/optical data. A survey of 7-ray pulsars in 
the IR can shed light on the ubiquity of the ultra- 
broad power-law magnetospheric component. 

Acknowledgements Based on HST programs 
GO-7836, GO-9156, GO-9797, GO-10600 and GO- 
11629). Support was provided by NASA through 
grants from the Space Telescope Science Insti- 
tute, which is operated by the Association of Uni- 
versities for Research in Astronomy, Inc., under 
NASA contract NAS5-26555. We also made use 
of archival Chandra observations provided by the 
HEASARC service. This work was partially sup- 
ported by NASA grant NNX09AC84G and by the 
Ministry of Education and Science of the Russian 
Federation (contract 11.G34.310001). We are in- 
debted to S. Zharikov for providing the VLT spec- 
trum. We also thank D. Thatte of STScI for pro- 
viding the dark frame required for the NICMOS 
reduction. 

REFERENCES 

Abdo, A. A., et al. 2010a, ApJ, 712, 1209 
— . 2010b, ApJS, 187, 460 

Anderson, S. B., Cordova, F. A., Pavlov, G. G., 
Robinson, C. R., & Thompson, Jr., R. J. 1993, 
ApJ, 414, 867 

Bignami, G. F., Caraveo, P. A., Mignani, R., Edel- 
stein, J., & Bowyer, S. 1996, ApJ, 456, L111+ 



9 



Brisken, W. F., Thorsctt, S. E., Golden, A., & 
Goss, W. M. 2003, ApJ, 593, L89 

Butler, R. F., Golden, A., & Shearer, A. 2002, 
A&A, 395, 845 

Caraveo, P. A., Bignami, G. F., & Mereghetti, S. 
1994, ApJ, 422, L87 

Chevalier, R. A. 1989, ApJ, 346, 847 

Chiaberge, M., Lim, P. L., Kozhurina-Platais, V., 

Sirianni, M., & Mack, J. 2009, Updated CTE 
photometric correction for WFC and HRC, 
Tech. rep. 

D'Angelo, C. R., & Spruit, H. C. 2011, MNRAS, 
416, 893 

Danilenko, A. A., Zyuzin, D. A., Shibanov, Y. A., 
& Zharikov, S. V. 2011, MNRAS, 415, 867 

De Luca, A., Caraveo, P. A., Mereghetti, S., Ne- 
groni, M., & Bignami, G. F. 2005, ApJ, 623, 
1051 

Greiveldinger, C, et al. 1996, ApJ, 465, L35+ 

Kargaltsev, O., & Pavlov, G. 2007, Ap&SS, 308, 
287 

Kargaltsev, O., Pavlov, G. G., & Romani, R. W. 
2004, ApJ, 602, 327 

Kargaltsev, O. Y., Pavlov, G. G., Zavlin, V. E., & 
Romani, R. W. 2005, ApJ, 625, 307 

Kern, B., Martin, C, Mazin, B., & Halpern, J. P. 
2003, ApJ, 597, 1049 

Koptsevich, A. B., Pavlov, G. G., Zharikov, S. V., 
Sokolov, V. v., Shibanov, Y. A., & Kurt, V. G. 
2001, A&A, 370, 1004 

Kuiper, L., Hermscn, W., Cusumano, G., Dichl, 
R., Schonfelder, V., Strong, A., Bennett, K., & 
McConnell, M. L. 2001, A&A, 378, 918 

Landsman, W., Bohlin, R. C, Neff, S. G., 
O'Connell, R. W., Roberts, M. S., Smith, 
A. M., & Stecher, T. P. 1998, AJ, 116, 789 

Marshall, H. L., & Schulz, N. S. 2002, ApJ, 574, 
377 

Mignani, R., Caraveo, P. A., & Bignami, G. F. 
1997, The Messenger, 87, 43 



Mignani, R. P. 2009, ArXiv e-prints, 0908.1010 

Mignani, R. P., Caraveo, P. A., & Bignami, G. F. 

1998, A&A, 332, L37 

Mignani, R. P., De Luca, A., & Caraveo, P. A. 
2000, ApJ, 543, 318 

Mignani, R. P., Zharikov, S., & Caraveo, P. A. 
2007, A&A, 473, 891 

Moon, D.-S., et al. 2004, ApJ, 610, L33 

Pavlov, G. G., Stringfellow, G. S., & Cordova, 
F. A. 1996, ApJ, 467, 370 

Pavlov, G. G., Welty, A. D., & Cordova, F. A. 
1997, ApJ, 489, L75+ 

Pavlov, G. G., Zavhn, V. E., & Sanwal, D. 2002, 
in Neutron Stars, Pulsars, and Supernova Rem- 
nants, ed. W. Becker, H. Lesch, & J. Triimper, 
273-+ 

Pavlov, G. G., Zavlin, V. E., Sanwal, D., Burwitz, 
v., & Garmire, G. P. 2001, ApJ, 552, L129 

Possenti, A., Mereghetti, S., & Colpi, M. 1996, 
A&A, 313, 565 

Rajagopal, M., & Romani, R. W. 1997, ApJ, 491, 
296 

Romani, R. W., Kargaltsev, O., & Pavlov, G. G. 
2005, ApJ, 627, 383 

Ruderman, M. 2003, ArXiv Astrophysics e-prints, 
10777 

Shibanov, Y. A., SoUerman, J., Lundqvist, P., 
Gull, T., & Lindler, D. 2005, A&A, 440, 693 

Shibanov, Y. A., et al. 2006, A&A, 448, 313 

The Fermi-LAT Collaboration. 2011, ArXiv e- 
prints 

Wang, F., Ruderman, M., Halpern, J. P., & Zhu, 
T. 1998, ApJ, 498, 373 

Wang, Z., Chakrabarty, D., & Kaplan, D. L. 2006, 

Nature, 440, 772 

Zavlin, V. E., & Pavlov, G. G. 2004, 
Mem. Soc. Astron. Italiana, 75, 458 

Zharikov, S., Mennickent, R. E., Shibanov, Y., & 
Komarova, V. 2007, Ap&SS, 308, 545 



10 



This 2-column preprint was prepared with the AAS lAlJ^X 
macros v5.2. 



11 



Table 3: 0656 FUV-MAMA total counts and fluxes in A-bins 



A-bin (A) Counts Background fx (10~^^erg/s/cni^/A) 



1153- 


1185 


48 


14 


177(37) 


1247- 


1282 


142 


33 


104(11) 


1313- 


1371 


190 


41 


88(8) 


1372- 


1497 


223 


46 


74(7) 


1498- 


1700 


114 


25 


57(8) 








Table 4 




COS FUV TOTAL COUNTS AND FLUXES OF B0656 


A (A) 


OA (A) 


Counts 


Background 


J" f-t n — 19 / / 2 / X \ 

JX (10 crg/s/cm /A) 


1130.1 


8.0 


170 


70 


138.8±18.5 


1138.0 


8.0 


228 


87 


176.6±19.4 


1146.0 


8.0 


231 


74 


177.5±17.9 


1153.9 


8.0 


214 


69 


149.0±15.5 


1161.9 


8.0 


244 


81 


154.2±15.0 


1169.8 


8.0 


234 


70 


142.8±13.6 


1177.8 


8.0 


189 


57 


107.2±11.3 


1185.8 


8.0 


223 


73 


114.2±11.4 


1193.7 


8.0 


241 


51 


137.5itll.6 


1247.1 


8.0 


273 


55 


129.0±9.9 


1255.0 


8.0 


274 


45 


132.3±9.8 


1263.0 


8.0 


203 


39 


92.7±8.2 


1270.9 


8.0 


222 


48 


97.5±8.4 


1278.9 


8.0 


251 


43 


116.6±9.1 


1286.9 


8.0 


260 


44 


122.0±9.3 


1319.5 


8.0 


195 


39 


94.6±8.5 


1327.5 


8.0 


144 


29 


71.2±7.5 


1335.5 


8.0 


205 


39 


105.5±9.2 


1343.4 


8.0 


198 


36 


105.3±9.3 


1351.4 


8.0 


187 


38 


98.8±9.1 


1359.4 


8.0 


184 


38 


99.4±9.3 


1367.4 


8.0 


182 


34 


103.0±9.5 


1375.3 


8.0 


120 


28 


65.0±7.7 


1383.3 


8.0 


171 


36 


98.2±9.5 


1391.3 


8.0 


121 


30 


67.6±8.1 


1399.2 


8.0 


161 


39 


92.7±9.6 


1407.2 


8.0 


147 


38 


84.4±9.3 


1415.2 


8.0 


146 


35 


88.1±9.5 


1423.2 


8.0 


121 


36 


69.2±8.7 


1431.1 


8.0 


136 


37 


82.4±9.5 


1445.0 


19.9 


308 


83 


78.3±6.0 


1464.9 


19.9 


249 


83 


61.5±5.6 


1484.9 


19.9 


286 


92 


77.5±6.5 


1504.9 


20.0 


236 


84 


65.6±6.3 


1524.8 


20.0 


210 


80 


60.8±6.3 


1544.8 


20.0 


234 


89 


73.6±7.3 


1566.1 


22.8 


232 


91 


69.0±6.9 


1593.4 


32.0 


293 


141 


.59.9±6.0 


1629.2 


40.0 


343 


173 


62.0±5.9 
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1800 2000 2200 2400 2600 2800 3000 

Wavelength, A 

Fig. 3. — The measured (absorbed) NUV-MAMA spectrum of B0656. The curve shows the best-fit absorbed 
PL model with ax = -3.01 for E{B -V) = 0.02. 



Table 5: B0656 NUV-MAMA total counts and fluxes in nine A-bins (+4 pixels shift in the dispersion direction 
is applied) 



A-bin (A) 


Counts 


Background 


fx (lO-i^erg/s/cmVA) 


1790-2040 


2594 


2316 


37(6) 


2040-2153 


1781 


1606 


27(6) 


2153-2295 


1889 


1634 


25(4) 


2295-2434 


1514 


1251 


22(3) 


2434-2608 


1590 


1274 


16.5(14) 


2608-2762 


1173 


909 


15.2(15) 


2762-2950 


1136 


913 


11.6(14) 
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Table 6 

acs/wfc ramp filter fluxes of b0656 



A (A) 


U (^Jy) 


4049 


0.31±0.3 


4401 


0.46±0.3 


4780 


0.38±0.3 


5165 


0.35±0.2 


5566 


0.45±0.3 


6244 


0.39±0.3 


6822 


0.49±0.2 


7393 


0.40±0.3 


7956 


0.48±0.3 


8751 


0.46±0.3 


9491 


0.35±0.5 



• HST broad 

• BTA 

• Subaru 

■ ACS ramp 



STIS 

F0RS2 

COS 
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Fig. 4. — Spectrum of PSR B0656+14 (as observed). The data plotted are a combination of fluxes in the 
literature (Table [2]) and our own work (Table [1]). 
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Fig. 5. — Spectral fits to the IR-UV spectrum of PSR B0656+14: continuum-only model (top), absorption 
lines (middle) or broad emission line (bottom) components. 
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Table 7 

Simple fits to the IR-UV spectrum 



Model 


Tcoid (K) 


civ 




fhot° 


(1014 Hz) 


X^/dof 


Continua 


500 


-0.29(6) 


0.404(8) 


0.37(2) 




86/40 


Absoption 




-0.32(5) 


0.424(9) 


0.37(2) 


2.91(3), 8.05(9) 


58/36 


Emission 


1200(400) 


0.4 


0.26(6) 


0.22(13) 


3.9(5) 


66/37 



Note. — Numbers in parentheses indicate uncertainty in the final digit (s). Where no un- 
certainty is given, the value is very poorly defined (< Icr). 

^Power-law {Fu = Npi^u~°'") normalization at u = 5 X 10^^ Hz, in ^tJy. 

''Flux density of the Rayleigh-Jeans power-law at v = 2 X 10^^ Hz, in fi3y. 



10^ 
IQi 

10" 

10-1 



fen 



IQ-l 
10-5 
IQ-'' 
10-"^ 




10" 1015 IQie IQi' IQi** IQi" 1020 1021 io23 io24 

MHz) 



Fig. 6. — Multiwavelength spectrum of B0656, including the NIR-UV data, Chandra X-ray fluxed spectrum 
(this work) and ^-ray fluxes from the 2FGL catalog. We also sketch the fit to the X-ray spectrum by 
Zavlin fc Pavlovl (|2004l) : two black-bodies with temperatures T = 0.82, 1.72 MK and a power-law with 
photon index T = 1.5. The X-ray absorbing column is A'h = 1 x 10^" cni^^, and the optical/UV reddening 
E{B -V)^ 0.03. 
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Table 8 

Pulsars with IR-UV and GeV detections 



Pulsar 


d 

(kpc) 


logr 


log Bo 


log Blc 


logE 


NIR 


Optical 


lo 

UV 


-^thcrm 




7-ray 


-Ecut 




Crab 


2.0(5) 


3.09 


12.58 


5.99 


38.7 


-5.3 


-4.8 


-4.2 




-3.9 


-3 


5.8(5) 




Vela 


0.287(18) 


4.05 


12.53 


4.65 


36.8 


-8.3 


-7.9 


-7.5 


-4.6 


-5.8 


-2.0 


3.2(1) 


1.30 


Geminga 


0.25(9) 


5.53 


12.21 


3.06 


34.5 


-7.0 


-6.7 


-5.7 


-5.3 


-4.3 


-0.1 


1.90(5) 


1.20 


B0656+14 


0.29(3) 


5.05 


12.67 


2.88 


34.6 


-6.6 


-6.3 


-5.6 


-2.8 


-4.8 


-2.0 


0.7(5) 


1.50 


B1055-52 


0.35(15) 


5.73 


12.03 


3.12 


34.0 




-6.3 


-5.4 


-3.4 


-4.5 


-0.85 


1.3(1) 


1.40 


B1951+32 


2.0(5) 


5.03 


11.69 


4.87 


36.6 




-5 




-4 


-3.9 


-1.7 


4.5(1.2) 


1.65 


J0437-4715 


0.156(1) 


9.83 


8.45 


4.45 


34.1 






-5.8 


-4.1 


-4.8 


-1.7 


1.3(7) 


1.55 



Note. — All values in typical standard units, age t in yr, magnetic field B in G, spin-down luminosity E in erg/s, and high-energy cut-off 
i?cut in GeV. The efficiency rj is the fraction of E seen in each spectral window. For the efficiencies, we evaluate vLi,/E at 2x10^'*, 5x10^**, 



2x10^^, and 1x10^* Hz, for IR, Optical, UV, and XpL respectively. For y^n^^j--^ we calculate the efficiency using th e bolometric luminosity of 
the hottest thermal component seen, and for the 7 efficiency we use the 100 MeV-100 GeV integrated luminosites in lAbd o et al. 1 l l2010bl '). The 
values of t, B and E have been corrected for the effect of proper motion, important for J0437— 4715. Fmw is PL photon index that connects 
the IR, X-rays and 7-rays, with typical uncertainty around 0.05 (Figure[7ll. 
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Fig. 7. — NIR/optical, X-ray and GeV spectra of pulsars detected in the NIR-UV window. Dotted lines are 
PLs drawn to connect the optical and 7-ray points (see text for discussion). Note that for J0437— 4715, the 
lowest frequencies are dominated by its WD companion (see Durant et al. 2011). 
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Fig. 8. — As Figure [71 but in terms of 77^ ee At:(PvFi,/E 
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A. Data reduction details 



A.l. STIS FUV grating 

The source was imaged on the Far- Ultraviolet Multi Anode Micro-channel Array (FUV-MAMA). The 
low-resolution grating G140L (which covers the wavelength interval « 1150-1700 A) with the 52" x 0'.'5 slit 
was used. The data were taken during two consecutive orbits (including the target acquisition). The total 
scientific exposure time was 4950 s|f|. FUV-MAMA was operated in TIME- TAG mode, which allows the 
photon arrival times to be recorded with 125 fj,s resolutiorlj. 

For each exposure, we processed the raw "high-resolution" images (2048 x 2048 pixels; plate scale of 0'.'0122 
per pixel — see §11 of the STIS IHB) using the most recent calibration files. As an output, we obtained 
flat-fielded low-resolution (1024 x 1024 pixels; plate scale 0'.'0244 pixel"^; spectral resolution 0.58 A pixel"^) 
images and used them for the spectral analysis. 

The processed images show a nonunif orm detector backgrou nd with a flat (constant) component and 
the so-called "thermal glow" component ( Landsman et al. 1998f) that dominates over most of the detector 
area and grows with increasing the temperature of the FUV-MAMA low-voltage power supply (LVPS) (the 
average LVPS temperatures were 42.60 and 42.84 C in the two consecutive orbits of our observation). The 
thermal glow is the strongest in the upper-left quadrant of the detector, where the dark count rate can exceed 
the nominal value, 6 x 10~^ counts s~^ pixel"^, by a factor of 80. To reduce the contamination caused by 
the thermal glow background, the source was imaged close to the bottom edge of the detector. 

We find the B0656's spectrum centered aX Y — 102 ± 2 pixels in each of the flat-fielded images (the 
centroid position slightly varies with X), where X and Y are the image coordinates along the dispersion and 
spatial axes, respectively. Even at this location on the detector the background still exceeds the nominal 
value by a factor of 1.5-3 (typical values are 1-2 xlO~^ counts s~^ pixel"^), depending on the position along 
the dispersion axis. To improve S/N, we co-added the exposures using the STSDAS task "mscombine" (the 
result is shown in Fig. [9]). The F-positions of the centroids differ by less than 3 pixels for different exposures 
and different wavelengths (V-positions). 

To subtract the strong, nonuniform background, we used a custom IDL routine with capabilities of 
grouping and fitting the backgroun d and selecting an optim al extraction box size depending on the position 
along the dispersion axis (see also iKargaltsev et al. 2004). The background is taken from the two strips, 
33 < y < 92 and 113 <Y < 172, adjacent to the source region, 93 < y < 112. To obtain the spectrum with 
a sufficiently high S/N, we have to bin the spectrum heavily; after some experimenting, we chose 5 spectral 



^The planned additional observations (6 orbits) were canceled because of the failure of STIS. 
^However, the absolute time of photon arrival is not known to better than 1 s (Shibanov et al. 2005). 







K lliTM— ! 











Fig. 9. — Raw FUV-MAMA spectrum of B0656. The boxes show approximate regions for the source and 
background extraction used in spectral analysis. 



20 



bins (A-bins; see Table [3|). The bins exclude the regions contaminated by the geocoronal emission (Lya line 
and the 01 line at 1304 A). Since the 01 line at 1356 A is essentially not seen, we do not exclude the region 
around 1356 A from the spectral analysis. The bins outside the contaminated regions were chosen to have 
comparable S/N (« 7-11), whenever possible. 

For each of the A-bins, we calculate the total number of counts, TVj, within the extraction boxes of different 
heights (one-dimensional apertures): As = 5, 7, 9, 7, and 7 pixels, centered at y = 102 for the first two A-bins 
and at y = 103 for the rest of the A-bins. To evaluate the background, we first clean the background strips 
(see above) from outstanding (> 10~'^ cts s~^ pixel"^) values ("bad pixels") by setting them to local average 
values. Then, for each of the A-bins, we fit the y-distribution of the background counts with a first-order 
polynomial (interpolating across the source region), estimate the number Nf, of background counts within 
the source extraction aperture As, and evaluate the number of source counts, Ns — Nt — Ni, (Table [3]). 

The uncertainty SNg of the source counts is evaluated as SNs = [Ns + SN^{1 + As/Ai,)]^/-^ , where SN^ is 
the Poisson background uncertainty scaled to the source aperture. We binned the distribution of background 
counts along the y-axis with the bin sizes equal to As and calculated SNh as the root-mean-square of the 
differences between the actual numbers of background counts in the bins and those obtained from the fit to 
the background. We calculated 6Ns and S/N for various extraction box heights and found the As values 
maximizing S/N for each A-bin (see Table [3]). 

We calculated the average spectral fluxes in the A-bins: 

(F ). -/aa, -^^^ ^■^^'^^ Ci ^^^^ 

where Ci is the source count rate in the i-th A-bin corrected for the finite size of the source extraction 
aperture, and R\ is the system response that includes the Optical Telescope Assembly throughput and 
accounts for the grating and slit losses and time-dependent sensitivity losses (Bohlin, 1999; see also §3.4.12 
of the HST Data Handbook for STI^ for details). 

A. 2. COS spectroscopy 

B0656 was observed with the Cosmic Origins Spectrograph (COS) aboard HST with the G140L grism for 
a total of approximately 22 ks. With the 1105 A central wavelength setting, the approximate sensitivity range 
is 1100-2250 A, with resolution R ss 5000 and cross-dispersion scale 0.023 "/pix. We used the standard, clear 
primary science aperture (PSA), a 2'.'5 diameter circle. The FUV-A detector was operated in TIME- TAG 
mode, so that the arrival time (with precision of 32 ms) and pulse-height were recorded for each photon, 
together with its position. The pulse height is not a measure of the photon energy, but can nevertheless be 
used for screening some spurious events by excluding values 1-4 and 3][1. The data were processed through 
the standard calcos pipeline, version 2.12, using the latest calibration files (January 2011) to produce 
corrected events lists (corrtag files). The raw 2-D image of the data is shown in Figure [101 Wavelength 
calibration spectra were taken during the observation, spatially offset from the source spectrum (Figure llOp . 

To extract fluxes from the data, we implemented an algorithm similar to the one used for STIS-FUV 
above. We converted corrected photon position to wavelength. Using the image in Figure llOi we found 
the source spectrum roughly along y — 495. Next, we split the data into bins of several columns each, we 
generated histograms of counts by y-coordinate in order to located the source centroid and measured the 
width by fitting a Gaussian plus constant. This provided an estimate of the position and of the extent of 
the spectral trace in the cross-dispersion direction as a function of wavelength, as well as a measure of the 
background. 

In each bin we summed the total number of counts falling within an extraction region, ±1.7 x ax pix about 
the trace, where a\ is the width of the fitted Gaussian profile. The value of 1.7 was chosen by experimentation 

^ |http: / /www. stsci.edu/hst / stis / documents/handbooks / currentDHB/STISJongdhbTOC.html 
' ^http : //www. stscl ■ edu/hst/cos/documents/lsrs/ 1 ISR2010_09(vl) .pdf 



21 



to maximize the signal-to-noise. This aperture captures 91% of the incident light. The measured flux in the 
aperture is then 

where N\ is the number of photons in the aperture, Bx is the fitted background surface density, and Ri^\ is 
the sensitivity function at position i in the dispersion direction within the bin. This flux was corrected to 
the standard calibration aperture (57pix wide), and multiplied by a factor of 1.03 to account for the slightly 
diminished sensitivity of COS at the time of observation compared to the benchmark of the calibration (see 
ISR COS 2010-ldl). We assume the uncertainty in the flux is dominated by Poissonian statistics in the 
source and background. Two regions of the spectrum could not be extracted due to strong geocoronal line 
contamination and the corresponding A-bins were omitted. 

A.3. STIS NUV prism 

We also retrieved the data for the program GO-9156 from the HST archive. This program consisted of 
two visits during which PSR B0656-I-14 was observed with STIS NUV-MAMA. The first and second visits 
occurred on 2001 September 1 (start time 52, 153.94 MJD) and 2001 November 16 (start time 52, 229.47 
MJD), respectively. The corresponding scientific exposure times were 6790.6 and 12761.0 seconds (3 and 5 
orbits, respectively). In these observations 52" x 0'.'5 slit together with the PRISM (with central wavelength 
of 2125 A, pixel scale 0'.'024) were used to obtain the dispersed spectrum of the target. NUV-MAMA was 
also operated in TIME- TAG mode. 

For the spectral analysis we used flat-fieldcd low-resolution images (1024 x 1024 pixels) and reference files 
provided by the STIS pipeline. We find the B0656's spectrum centered &iY — 499 ± 2 pixels in each of the 
fiat- fielded images (the centroid position slightly varies along the dispersion axis X). To increase S/N we 
combined 5 exposures from the second (longer) visit using STSDAS mscombine task. The resulting image is 
shown in Figure [TTJ Most of the emission is confined between X = 490 and X = 570. 

To extract the NUV spectrum, we use the procedure similar to that applied to the FUV data. We scan 
the count distribution within two strips, 430 < 1" < 489 and 510 <Y < 569, adjacent to the source region, 
490 <Y < 509. The average background count rate at the source location is 1.2 x 10~^ counts s^^ pixel"^, 
i.e., about two orders of magnitude higher than in the case of MAMA-FUV observation. Therefore, to obtain 
the spectrum with a sufficiently high S/N, we bin the spectrum into 8 spectral bins with S/Nrs 5-11. 

Contrary to grating observations, the conversion from the pixel number (along the dispersion axis) to 
wavelength is nonlinear for PRISM spectrgF^I. Since the PRISM dispersion depends strongly on the wave- 
length (the resolution decreases with increasing wavelengtt{3) , even a small shift (in terms of pixels) will 
translate into a large shift in terms of the assigned wavelength at longer wavelengths (A > 3000 A). Fur- 
thermore, because of the strong dependence of the PRISM throughput on wavelength, assigning a slightly 
incorrect wavelength can result in a large error in the derived fiux at A > 3000 A. Hence, it is important 
to accurately determine the wavelength zeropoint. The latter depends on the actual position with respect 
to the aperture (slit) center as well as on the offsets caused by the Mode Selection Mechanism (MSM) 
(positioning/tilt of the slit/PRISM) non- repeatability and thermal drift. 

The wavecal exposures taken between the science exposures are expected to accurately account for the last 
two effects while the source position within the slit can be measured from the CCD images taken during each 
visit through the same 52" x 0.5" slit aperture. Although at the first glance the source appears to be well 
centered within the aperture, a more detailed investigation revealed a 0"045±0"02 shift (about 9% of the slit 



^http : //www ■ stsc l ■ edu/hst/cos/document s/ | isrs/ISR2010_15(vl) .pdf 

l^See §12.1 of the STIS Instrument Handbook for details 
"See §4.4 of the STIS IHB ( http : //www . st s c i . edu/ , 

hst/st is/document s/handbooks/currentlHB/ 

c04_spectros5 . html#310764) 
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width) of the position of the best-fit source centroid with respect to the sUt center. This relatively small shift 
amounts to 1.6 MAMA-NUV pixels. Such a difference in assigning the MAMA-NUV wavelength zeropoint 
would produce a noticeable difference in the flux at A > 2600 A because the throughput curve is steep and 
the resolution very low for these wavelengths. Hence we account for this shift when assigning the wavelength 
zeropoint. Since the measured shift has a large associated uncertainty, we choose to ignore the spectrum 
at A > 2950 A (at these wavelengths the 0'.'02 uncertainty may have a large impact of the flux calibration). 
Additionally, the steep throughput curve and low dispersion (compared with the line-spread function) mean 
that the apparent flux at wavelengths A >3000 A would be contaminated by shorter wavelengths. This effect 
is difficult to account for, giving another reason to consider the flux calibration unreliable for A > 2950 A 
and limit the extraction of the NUV spectrum to the 1790 — 2950 A range. 

Aside from the caveats above, we use the same flux calibration procedure as for the STIS FUV data to 
derive the NUV spectrum. 

A.4. ACS/WFC Ramp Filter Photometry 

The Advanced Camera for Surveys (ACS) aboard HST provides three ramp filters for narrow-band imag- 
ing. A ramp filter has a band-pass which depends on the part of the filter that light passes through, hence 
one can choose from a range of central wavelengths, with about 10% bandwidth in each case. We utilized 
each of the three ramp filters, for a total of 11 photometric points, see Table [TJ The pipeline-produced, 
drizzle-combined images are shown in Figure 1121 

We use the same aperture size for all filters and central wavelength settings, because the width of the 
PSF varies only slightly with wavelength (^^15% across the optical range). We used r = 0'.'125 (2.5 pixels), 
which is close to the typical FWHM, and gives nearly optimal SNR in all the images. The encircled energy 
fraction within the aperture varies from 52% in the red to 70% in the blue. The background was taken from 
an annulus with the same center and radius 1"< r <1'.'75. To estimate the aperture correction, we measured 
the fiux within r — 0'.'125 and 0'.'5 (the stand ard calibration apertur e), for the nearby brighter star N7, 



tne nux witnm r — u: izo ana U. o [the stand ard can oration apertur e), tor tne nearoy Drignter star 
located approximately 10" northwest of B0656 ( Koptsevich et al1l200l[ ). This is one of the few field sources 



that has a point-source profile (most field sources are extended, i.e., galaxies), and photometry shows that it 
is indeed a star (see Section lA^ . Since the filter band-pass is position-dependent, it is necessary to choose 
a calibration source for the aperture correction as close to the source of interest as possible. Finally, the 
fluxes were calculated using the photometric zero point delivered by the pipelinj^. 



The CTE correction was implemented following the prescription of IChiaberge et al.l (|2009lj^4 We in- 
terpolated between the August 2005 and March 2006 values to find the parameters appropriate to our 
epoch of observation, giving a mean correction ^mag — x SKY'' x FLUX'^ x (Ytran/2000)=0.04mag, 
where a = 0.40,6 = — 0.26,c = —0.47, and the number of transfers in the y direction during readout was 
Ytran «1100 for 4049, 4401, 4780, 5165, 7956, 8751, 9491 A and Ytra„ -990 for 5566, 6244, 6822, 7393 A. 

A.5. NICMOS photometry 

We downloaded the NICMOS data and the best available calibration files from the archive. Each image is 
a 636.1 s exposure in ACCUM mode. This mode is non-standard: there are multiple reads at the end of the 
exposure rather than reads at regular intervals throughout the exposure. This leads to inefficient screening 
of CRs. Four images per orbit were obtained using a dither pattern, with one orbit for FllOW (2554s total 
integrations), two for F160W (5089s total integration) and three for F187W (7633s total integration - see 
Table [ID . 

We processed the data, utilizing the custom-made dark file provided by STScI at our request. This was 
required due to the time-variability of the detector temperature both with the age of the instrument and 



^^www. stsci . edu/hst/acs/documents/isrs/isrOTll .pdf 
^• ^http: //www, stscl ■ edu/hst/acs/dociimeiits/isrs/lsr0901 .pdf | 
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time ellapsed since its previous power-on. The processing of the raw data into final images shown in Figure 
[13] included the following steps: 



• the pedestal correction removes bias-like additive signal from each amplifier quadrant by calculating 
the median pixel values, after rejection of stars and hot pixels; 

• the task multidrizzle resamples images to a common reference frame, removing geometric distortion. 
It also compares the median of these images with the input resampled images, in order to identify and 
reject some cosmic rays/bad pixels; 

• final images were made by a median combination. We found that the median gave better result than 
the average, since the ACCUM mode results in a very high incidence of cosmic rays. 

Photometry was performed using small apertures of r = 2.5 pix (0'.'19), with the aperture correction to 
the O'.'S nominal aperture calculated from the photometry of Star N7. The dither pattern sometimes took N7 
close to the chip edge, so we used only those frames in which N7 was sufficiently far from the edge to avoid 
low sensitivity pixels in the O'.'S aperture. The calibrated zero point is defined relative to 1.15 times the flux 
in the 0'.'5 aperture (equivalent to correcting to the "infinite aperture"). The final flux was calculated using 
the header PHOTFNU keyword, provided by the pipeline. 



A. 6. Optical-IR calibration check 

In an effort to check the calibration of our observations, we consider the flux measurements of field star 
NT (see Figure [H]), also known as 2MASS J06594760-I-1414253, in the nominal-size aperture 0'.'5. Using 
broad-band photometryfiJST in the optical and 2MASS in the NIR) to estimate the spectral class and 
reddening of the sourcqlj, we find that the star spectrum approximately matches a 4000 K (K-type) giant 
at 2.5 kpc (i.e., the Perseus Arm) reddened by E{B — V) — 0.35 (the photometry also matches a dwarf at 
~80pc, but the reddening is unreasonably high for such a small distance). The broad-band fluxes and the 
fluxes we derive appear to match very well in the optical. 



Our IR fluxes for N7 are about 10% higher than those of iKoptsevich et al.l (j200ll ). The likely reason for 
the discrepancy is that that these authors also used the exposures of N7 which fell onto the low-sensitivity 
edges of the detector. Our IR fluxes are, however, about 20% lower than the 2MASS fluxes. The 2MASS 
extraction aperture includes the fluxes of three faint nearby sources (resolved by ffST), which contribute 
roughly 10% to the NIR flux each, and give summed fluxes consistent with the 2MASS ones, within the 
uncertainties. 



See our star type fitting tool at ]http : //www ■ astro ■ uf 1 ■ edu/] 
~martin. durant/Sclass .html 
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Fig. 10. — COS image of PSR B0656+14, showing events filtered by pulse height and data quality (binning 
for display purposes is 5x2 original pixels, 0'.'115x0'.'184). Approximate spectral and sky extraction regions 
are shown (see text). The bright strip across the top of the image is the calibration spectrum, and the two 
prominent bright patches are images of the round aperture centered on bright geocoronal emission lines, 
which look elliptical after binning. 
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Fig. 11. — Raw NUV-MAMA spectrum of B0656. The boxes show approximate regions for the source 
background extraction used in spectral analysis. 
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Fig. 14. — Optical-IR spectrum of field star N7, showing 2MASS and HST broad-band photometry (black 
points), NICMOS photometry (red squares) and ACS narrow-band photometry (blue circles). Statistical 
uncertainties are quite smal in most bins, and hence not seen. 
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